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ABSTRACT 

The mode of explosive burning in Type la SNe remains an outstanding problem. It is generally thought to 
begin as a subsonic deflagration, but this may transition into a supersonic detonation (the DDT). We argue that 
this transition leads to a breakout shock, which would provide the first unambiguous evidence that DDTs occur. 
Its main features are a hard X-ray flash (~ 20 keV) lasting ~ 10~ 2 s with a total radiated energy of ~ 10 40 ergs, 
followed by a cooling tail. This creates a distinct feature in the visual light curve, which is separate from the 
nickel decay. This cooling tail has a maximum absolute visual magnitude of M v as —9 to —10 at » 1 day, which 
depends most sensitively on the white dwarf radius at the time of the DDT. As the thermal diffusion wave 
moves in, the composition of these surface layers may be imprinted as spectral features, which would help 
to discern between SN la progenitor models. Since this feature should accompany every SNe la, future deep 
surveys (e.g., m = 24) will see it out to a distance of « 80 Mpc, giving a maximum rate of ~ 60 yr _I . Archival 
data sets can also be used to study the early rise dictated by the shock heating (at « 20 days before maximum 
B-band light). A similar and slightly brighter event may also accompany core bounce during the accretion 
induced collapse to a neutron star, but with a lower occurrence rate. 

Subject headings: hydrodynamics — shock waves — supernovae: general — white dwarfs 



1. INTRODUCTION 

The use of Type la supernovae (SNe la) as cosmological 
distance indicators has brought attention to the theoretical 
uncertainties that remain about these events. Generally it is 
thought that they result from the unstable thermonuclear igni- 
tion of a C/O white dwarf (WD). One key question is how 
the burning front propagates during the incineration. The 
consensus is that the flam e begins as a subsonic deflagration 
dNomoto et al.ll 197611 19841) to match the observed nucleosyn- 
thesis and light curves (Filippenko 1997), but the later propa- 
gation is more uncertain. Motivated by terrestrial combustion, 
many have argued for a delayed detonatio n transition (DDT, 
lKhokhlovll99UIWooslev & Weaveri[l994l) . Alth ough a DDT 
may be needed to bette r replicate observations dPlewa et al.l 
120041: iLivne et all 12005b. how and if it happ ens is uncertain 
dNiemever & Wooslevll 9971; iNiemeverl 1 999; Wo oslevl2007b . 

Any observational feature that unambiguously demon- 
strates that a detonation occurs would be helpful for resolving 
this uncertainty. One consequence of a detonation is that it 
drives a shock through the WD surface layers, the breakout 
of which is a clear signature of the DDT. For Type II super- 
novae (SNe II), it was long expected that shock breakout from 
core-bounce would produce an X-ray and/or ultra-vi olet flash 
(IColgatdll974l IKlein & Chevaiiel[l978T: lFal^[T978h . As the 
earliest electromagnetic emission available, it is an important 
probe of the progenitor st ar and its circumstellar environment 
(Matzner & McKee 1999). These X -ra y flashes have been ob - 
served in the cases of GRB 060218 (Cam pana et al.|p2 006). 
which was associated with the SN Ic 2006aj (Mazzali et al. 
120061) and t he X-ray transient XRT 080109, associated with 
SN 2008D dSoderberg et al.l 12008}) . Indeed, since the SN 
la rate in star formi ng regions is ~ 1/3 that of SNe II 
dMannucci et al.ll2005b . there may be hope that a SN la shock 
breakout may soon be seen. The question is how bright is 
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such an event and how long does it last 

In this paper, we characterize the light curve expected from 
an SN la shock breakout. In $2j we present the initial WD 
profile before shock passage and investigate when the shock 
runs away form the detonation using both a numerical, hydro- 
dynamic simulation and analytic arguments. We then present 
equations that describe how the shock steepens as it propa- 
gates through the surface layers ahead of the detonation. In 
$3] we follow the expansion and cooling of the shock-headed 
envelope using a semi-analytic analysis. We demonstrate that 
the entire light curve can be modeled with self-similar solu- 
tions. The breakout produces a sharp, hard X-ray flash with 
a duration of w 10~ 2 s (also see Imshennik et al.|[198 11) . As 
the shocked envelope expands and cools, the lightcurve shifts 
from the X-ray to the UV and eventually to the visual wave- 
bands. In $4] we show that it reaches a peak absolute mag- 
nitude of My = -9 to -10 at ss 1 day. We also discuss what 
depths of the WD are probed by the receding thermal diffu- 
sion wave as a function of time. In $5] we summarize our 
results and discuss the possibility of shock breakouts from ac- 
cretion induced collapse (AIC) to a neutron star. 

2. THE DYNAMICS AND THERMODYNAMIC 
PROPERTIES OF THE SHOCKED ATMOSPHERE 

We first summarize the properties of the unshocked WD 
envelope ( jj2.1| i and then model the detonation to determine 
how the outward propagating sho ck st eepens as it runs down 
the background density gradient ( §2.21 . From these we derive 
the equations that desc ribe how the envelope responds to the 
shock's passage Q2.3b . 

2.1. The Pre-shocked Atmosphere 

The pre-shocked envelope, which we denote with the sub- 
script 0, is estimated as a constant flux atmosphere, F = 
(TsBTrff' an d i s m radiative equilibrium. We assume a Chan- 
drasekhar mass (M = 1.4M Q ) WD with a radius « (3 - 
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6) x 10 8 cm, which is larger than the pre-ignition radius due 
to expansion during the deflagration phase (in contrast, the 
pre-ignition radius is « 1.6 x 10 8 cm). Radiative diffusion 
gives 



1 eff- 



3k dPu 



(1) 



where we take a constant opacity, k = 0.2 cm 2 g corre- 
sponding to electron scattering in hydrogen deficient mate- 
rial. We take the gravitational acceleration g = GMjR\ to be 
constant in these shallow, surface layers. Integrating equation 
{T), assuming an ideal gas equation of state, Pq = Po^b To / pm p , 
where fi is the mean molecular weight, we get the pressure as 
a function of density, 



Po( A) ) = 6.1xl0 1 



-1/3-4/3 4/3 _3 



(2) 



where gg = g/10 9 cm s" 2 , = r e ff/10 5 K, po is in cgs units, 
and [1 = 4/3 (for a helium dominated composition). 

At sufficiently large depths, the pressure becomes domi- 
nated by non-relativistic electrons, in which case 



Po(po) = 9.91 x 10 I2 (p /^) 5/3 ergs cm" 



(3) 



where fi e is the mean molecular weight per electron. This 
switch occurs around po w 9 x 10 3 g cm -3 . These two power 
laws (eqs. J2) and Q) motivate our use of a polytropic back- 
ground model in the following sections. 

2.2. Detonation and Shock Runaway 

The energy budget available to the shock is set by where it 
runs ahead of the detonation, giving a characteristic velocity 
and density for the shock, V lun and p mn . The strong shock sub- 
sequently steepens in the declining density gradient according 
to 



V = V m 



Po 

pmn 



-0 



(4) 



with = 0.18 58 for a radiation pressure dominated shock 
(lSakurailll96a -/ = 4/3). 

We use a one-dimensional, hydrodynamics code (described 
previously in Weinberg & Bildsten 2007) to follow the det- 
onation's propagation and transition to a shock, which we 
show in Figure Q] The background is initialized with a con- 
stant flux radiative atmosphere, which smoothly transitions 
to a nearly isothermal degenerate star. A detonation is initi- 
atedat 10 7 g cm" 3 (near where the DDT is expected, Wooslev 
2007) by artificially raising the temperature. Above this depth 
the composition is a flammable mixture of equal parts 12 C 
and 16 0, and below it is inert material to prevent the deto- 
nation from propagating inward (in accord with the preced- 
ing deflagration). The simulation shows an outward prop- 
agating detonation, which is accompanied by a weak shock 
sent downward into the WD core. The detonation propagates 
until w 2 x 10 6 g cm" 3 , at which point the shock runs away. 
This is most clearly seen in the top panel of Figure Q] where 
the synthesis of intermediate mass elements falls off in com- 
parison to la O (the depletion of 12 C at these depths is due to 
residual burning behind the shock). In the bottom panel, we 
see that the velocity falls off as a Vo oc /Oq 0186 power law as 
expected (eq. HJ). From this plot we determine that the char- 
acteristic numbers at runaway are p mn w 2 x 10 6 g cm" 3 and 




Initial Density (g cm ) 

FIG. 1 . — The composition (top panel) and fluid velocity (bottom panel) as 
a function of initial WD density. The detonation is initiated at 10 7 g cnr 3 , 
sending a burning wave into the flammable lower-density region, and a shock 
wave down into the WD. At 2 X 10 6 g cm 4 , the steepening shock begins 
to race ahead of the burning, as can be seen by the marked decline in in- 
termediate mass element synthesis. At this depth, the shock has a speed of 
6 X 10 s cm s" , which sets the initial conditions for the shock breakout 
calculation. 



Similar results are obtained by using simple analytic scal- 
ings. As the detonation propagates, it causes a change of pres- 
sure 



/?det 

Po 



(5) 



where E nuc rs 0.8 MeV nucl 1 is the specific energy release 
from burning two 12 C and two ls O to 56 Ni, E mt = 3Ep/4p e is 
the internal energy, and 



E F = 0.41pg /3 (2/^,) 1/3 MeV nucl" 



(6) 



is the Fermi energy. The speed and temperature of the det- 
onation are then r det = (3p d et/a) 1/4 and V det = (pda/ Po) 1/2 c s , 
respectively, where c s is the sound speed. At the depths where 
the electrons are degenerate and relativistic, this gives a den- 
sity independent detonation speed of Vdet ~ 6 x 10 8 cm s" 1 , 
consistent with our numerical results. 

The burning in the detonation occurs over an "induction 
length" given by 



A = (£nuc/e)V df 



(7) 



Vn 



i 6 x 10 cm s 



where e is the energy generation rate. As the density de- 
creases, the burning rates does as well, increasing A. Once 
A is greater than a pressure scaleheight H = P/ pg, the burning 
can no longer keep up and the shock races ahead of the burn- 
ing front. Using the energy generation rate, e, for burning 12 C 
to 24 Mg from Woosley at al. (2004, see their eq. [4]), this 
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occurs at a critical density of 

p run «2xl0 6 ^ 11 gcm- 3 . 



(8) 



The prefactor is robust because of the strong temperature and 
density scalings of the energy generation rate. This compares 
well with our numerical calculation (see Fig. [TJ. 

2.3. Shock Heating and Adiabatic Expansion 

The jump conditions for a strong shock (Whitham 1958, 
his eq.[32]) give the pressure p^, density p sn , and velocity V 
behind the shock in terms of the shock velocity Vq, 



Psh : 



7+1 



PoV , p s h : 



7+ 1 

— T 



7- 



7+ 1 



V . 



(9) 



We next need to understand how these properties change with 
expansion. 

Consider a shell of mass AM, density po, and thickness 
Aro <C /*o, where ro w is the local radius in the WD enve- 
lope. The shell's radius as a function of time is 



r(p ,t) = r + V(p )t « R* + V(p )t . 



(10) 



Pressure gradient s may increase the veloc ity V (po) by as much 
as a factor of 2 (Matzner & McKee 1999), but we find this 
difference has little effect on our solutions. By continuity of 
mass, the shock compresses the shell to a thickness Ar sn = 
(7- l)/(7+ l)Aro. Subsequently, the thickness increases due 
to the difference in velocity between the top and bottom of the 
shell, giving a thickness as a function of time 



Ar(p ,t) = Ar sh + AV(p )t, 



where AV{po) is 



AV(po) 1 



dV 
dpo 



dp 



Ar . 



(11) 



(12) 



We take the equation of hydrostatic balance, dPo/dro = -pog, 

and assume a polytropic equation of state, Po = Xpo" + 
(where K and n are set by either eq. 12 or 0). Setting 

Aro ~ Ho = Po/ pog = Kp\j n /g, where Hq is the pressure scale 
height, we find 



AV(p Q ) = 



n/3 



n+ 1 7+ 1 



-V™ 



Po 

Prun 



(13) 



which is smaller than V by a constant factor of n/3/(n+ 1) : 
0.14 (forn = 3 and /3 = 0.186). 
By mass conservation 



4np Rl Ar = 4tt p(p , t) r(p ,t) 2 Ar(po,t), 



(14) 



where p(po,t) is the expanded density as a function of time. 
This gives us a relation for the fractional change in density 



p(po,0 
Po 



Aro 



r(p Q ,t) 2 Ar(po,ty 



(15) 



which is a function of only po and t . 

From the shock jump conditions (eq.j9)), we find the pres- 
sure in the shocked and expanded shell, 



Psh(po) = —tPoVq = — — - 
7+ 1 7+ 1 



1-2/3 



— • (16) 



Po 

Prun 



We assume that the pressure decreases adiabatically as the 
shell expands, i.e., P oc p 1 , so that 



P(po,0 = Psh(po) 
2 



p(Po,0 



Po 



7+ 1 



PrunKn,, 



(—) 


p(Po,0" 


\ Prun J 


Po 



1 T 



(17) 



Although the adiabatic approximation is goo d, it is in general 
not exactly true, an issue that we discuss in 33.31 By combin- 
ing equations dl3T > and ( [T7l i. we solve for the pressure of the 
expanded shell with initial density po at any time t. 

3. THE LIGHT CURVE FROM SHOCK BREAKOUT 
AND SHOCK-HEATED COOLING 

3.1. The Prompt Breakout Flash 

Photons in the radiative shock stream out once r < c/V, 
where r is the optical depth. Substituting V —> TV, where T = 
[1 -(V/c) 2 ]~ 1//2 (the shock is mildly relativistic at breakout), 
we find a breakout depth of 



P0,br = Pn 



7+1 Tc g 



2 V™ 



KKpt 



l+l/n 



l/(l+l/«-/3) 



r^m-S/' Iv \0.87x/-0.87 -0.16ap /1 on0.87 „ -3 

~10 (g9/Ai3.8) V 9 p 6 (Iy 1.7) gem , 

(18) 

where ^i 3 . 8 = K/6 x 10 13 cgs, V 9 = V mn /10 9 cm s~\ and p 6 = 
Prun/10 6 g cm" 3 (using typical values from Fig.|TJ. 

This initial streaming of photons gives rise to the prompt 
breakout flash. The energy density available is E = 3p s h, and 
using equation ( fTBI ) for /? sn , we estimate a total energy budget 
at a depth po,br of 



£flash = 4ttRIH x 3p s 



sh 



'A x 10 40 (gg/Kn.sT^V^p^Rts ergs, (19) 
where R% 5 =R*/3 x 10 8 cm, and associated temperature 

7flash = (3p sh /fl) 1/4 

«2xl0 8 fe/^3. 8 ) 01 V 9 032 Pr 68 K. (20) 

This temperature is a lower limit, since non-LTE effects 
(which we ig nore to simplify our analysis) will only harden 
the spectrum dKatz et al.l |2 009). We further discuss the limita- 
tions of this assumption in 33. 6| These numb ers are consistent 
with the findings of Imshenni k et alj dl981l) . who predicted a 
hard (m 20 keV) X-ray flash associated with WD detonation. 
The main difference is that we are quantifying this flash in 
terms of the DDT picture. The timescale for the energy re- 
lease will be dominated light travel effe cts ( R*/c w 10~ 2 s), 
which we include in our calculations in 33.21 Also, the DDT 
will not occur ever ywhere at once (fo r example, as found in 
the simulations by Rtipke et al. 2007), which further smears 
out the breakout flash. 

3.2. Cooling of the Ejected Shock-Heated Envelope 

Following breakout, a thermal diffusion wave begins prop- 
agating back into the shock-heated envelope, releasing pho- 
tons that will be seen as the shock breakout flash. At any 
given depth in the shock-heated envelope, the energy density 
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Igio t ( s ) 

FIG. 2. — Solutions for the cooling, shock-heated WD envelope with V rlln = 
10' cm s - ' and p run = 10 6 g cm -3 . The top and middle panels show the time- 
dependent luminosity and effective temperature from the moment of shock 
breakout until rs 1 day. The bottom panel shows the thermal diffusion depth 
Arjiff = Ar(t = fdiff). The solid lines (dashed lines) are for an initial WD 
radius of R* = 3 X 10 s cm (6 X 10 s cm). In the bottom panel, the horizontal 
dotted line denotes i?* = 3 X 10 s cm. 

is E = 3p and leaks out of the envelope on a thermal diffusion 
timescale, which at a depth po and a time t, is 



fdiff(/0(),f)= ~p(po,t)[Ar(p ,t)] 
c 



(21) 



Since t$s sets the timescale when an observer sees down to a 
certain depth in the shock heated envelope, we set fdiff (po,f) = 
/ to solve for po(t), using the prescription for p(po,t) and 
Ar(po,t) from 32.31 Once we have po(t), we can solve for 
any other property of the envelope as a function of time. In 
particular, the luminosity of the expanding and cooling enve- 
lope is 



my. 



Airr{t) 2 E{t)c 



where 



r{t) = nAr(t)p(t) 



(22) 



(23) 



is the optical depth of the diffusion layer. 

The general solution of Lit) is non-algebraic, but easily 
solved numerically. In Figure [2] we plot L(t) (top panel), the 
effective temperature {middle panel) 



T e «(t): 



Lit) 



"I 1/4 



47rr(f) 2 crsB 



(24) 



and the thermal diffusion depth Ar^ = Ar(f = f^ff) {bottom 
panel). The luminosity decreases as a broken power law in 
time. The initial decline is steeper, with a L oc f~ 4 / 3 power 
law. In practice, this power law and the very high initial lu- 
minosities are an artifact of our one-dimensional treatment, 
and will not be seen in observations. Just as for the breakout 
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FIG. 3. — The early time light curve, including time travel effects across 
the WD surface (thick lines) in comparison to the one-dimensional solutions 
plotted in Figuref2](</iM lines). Light travel effects decrease the observed lu- 
minosity at early times when light from the entire surface has not yet reached 
the observer. 



flash, light travel effects and not simultaneous DDT ignition 
will smear out the light curve over R*/c w 10~ 2 s or more. To 
better quantify this effect, we plot the light curves altered by 
light travel effects in Figure[3] These are calculated according 
to 



Ut)- 



•rr/2 



L[t r {t , cos 9)]d cos 6, 



(25) 



where t r = t- r{t r ){l — cos 9) jc is the retarded time of photons 
emitted at a angle with respect to the observer. Note that 
r{t r ) must be evaluated at the retarded time because the ra- 
dius was smaller in the past The robust feature to take away is 
that - 10 40 ergs of energy is released during this initial peak 
(consistent with eq. fl9l ). In the second stage, the luminos- 
ity decreases more slowly, from 10 42 ergs s" 1 to 10 40 ergs s -1 
over ~ 1 day. The effective temperature evolution covers a 
wide range of wavelengths, going from X-rays to ultraviolet 
to visual over the course of w 1 day. 

The power law behavior of L{t) can be understood ana- 
lytically by considering the expansion and cooling in cer- 
tain limiting cases. To emphasize this, we have plotted = 
3 x 10 8 cm in the bottom panel of Figure [2] {dotted line). A 
comparison with L(t) clearly shows that the luminosity break 
occurs when Ar^f has expanded to a thickness ks R*. Moti- 
vated by this, we consider the time-dependent behavior in two 
different limits: 

1. At early times, t <C R*/V, the envelope has barely 
moved and r{t) w R*. By equation ( fT5t , we see that 
p/po(xt~ l since the radial position of the shell remains 
constant and the thickness increases as Ar oc t. 

2. At late times, t ^ R*/V, both r and Ar increase like f, 
so that pj po oc f~ 3 . 

We next derive the analytic self-similar solutions in each limit. 
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3.3. Case 1: t <.R*/V 

In this limit the expansion is plane-parallel so that the ra- 
dius is basically fixed at r(t) w R*, whereas the thickness has 
expanded by Ar(t) i=s AVt . Since p(t) oc t , we find that the 
thermal diffusion depth is independent of time when we set 
the condition t^s = t, and we find 



PO.diff : 



AV 



;4x 10" 



(g9/£i3.8) v 9 p 6 gem . (26) 

The physics expressed by equation d26| > is that the thermal dif- 
fusion wave sits at a fixed depth because the column of mate- 
rial sitting above it is fixed. The optical depth is proportional 
to column, and thus also fixed, 



r diff = nAVtp(t) = nHopo rj l{g 9 /K n f*% 

This gives a luminosity of 
T 4TrRlE(t)c 



0.87 T /- 0.87 -0.16 



P6 



(27) 



<9xl0« to/^i3.8)-°- 34 V 9 °- 72 p° 13 ^ 



,-4/3 

f_ 2 ergs s 



(28) 



where f_2 = t / 10 2 s and 

r eff « l x io 7 (g 9 /«i 3 .g 



r 0.086 y 0.1 8jO 0.033 i: l/3 K (29) 



is the effective temperature. 

Since in this limit, the same mass shell is always doing the 
radiating (as indicated by a fixed po.diff m ec l- EH), it can- 
not be evolving adiabatically (thus our assumption in eq. fTTl 
is not strictly valid). In a more detailed study we solve the 
entropy equation in the plane parallel limit, including the ra- 
diative loss term (Chang, Matzner, & Piro, in preparation). 
We find that the cooling is only marginally faster when we 
include radiative losses, with a logarithmic time dependence. 

The plane-parallel limit applies until Vt w R*. Setting 
Po.diff = Po mto equations and ||9), we find this occurs at 



\0.16 T/ -1.2 -0.22 



^8.5 S 



(30) 



t « 1 x 10" 2 ( 8 9/K l3 .,f lb V 9 --p b ' 

at which point we must begin to account for the radial expan- 
sion. 

3.4. Case 2: r »/?»/V 

The layer has now moved an appre ciable distance fr om the 
star, so that r(t) » Vt (as addressed in lChevaiierl lT992 ) . Solv- 
ing for the diffusion depth in this limit 



PO.diff = Pn 



n + 1 2 



Vm 



n P 7+1 K Kp, 



l+l/n 
run 



~2 (g 9 /Ki 3 .&) 
The optical depth is 



0.66t/0.66 0.12 d -1.3.1.3 



R 



8.5 



Ct 

R* 

gem" 



l/U+ZJ+l/n) 



(31) 



c 
AV 



n+1 7+ 1 c 

nfi 2 V mn 



T, / /p- s.0. 12, 7-0.88 
= 21(^9/^13.8) V 9 p 6 



Po 

Prun 
0.16 



D-0.25,0.25 



This gives a luminosity of 



L{t) « 3 x 10- ig 9 /K n ^V^p^R^ A ergs s 



(32) 



(33) 



which has a power law time dependence, L oc t~° M , consistent 
with what IChevalierl dl992h found in the same regime. The 
effective temperature is 

7k(o « i x io 6 (g 9 /^3.8)- 0065 v 9 001 vr 035 ^:^^" 046 k. 

(34) 

Since the surface area is increasing and the diffusion wave 
is moving into higher temperature material, the evolution is 
shallower than in the plane-parallel case. 

3.5. The Exposure of Originally Degenerate Material 

At sufficiently late times, the diffusion wave moves into 
material that was originally degenerate before being hit by 
the shock. Setting po.diff in equation OTb equal to w 9 x 
10 3 g cm" 3 ( 92. 11 ). this occurs at w 700 s. The dynamics are 
essentially the same as that of §3.2, but now n = 3/2 instead of 
n = 3. The diffusion wave now moves through the layer more 
slowly, with pdiff.o ocf , which gives a luminosity, 



r W ~1v in 40 <„ IV \-0.41 v 1.9 n.36p0.83,-0.16 \ 

L(t)ttZX W (g 9 /Ki3) V 9 p 6 K S5 ? 4 eigss , 



(35) 



and effective temperature, 

T M-lvin 4 r„ /if \-0.05 8 T70.030 f).0058p0. 11,-0.44 v 
hmt)- 2x 10 {g 9 /Ki3) V 9 p 6 K S 5 t 4 K, 

(36) 

where = K/10 13 cgs and f 4 = f/10 4 s. The initially de- 
generate material causes the light curve to flatten and the lu- 
minosity remains at ~ 10 40 ergs s" 1 well until a day after the 
initial shock breakout (as is shown in Fig. |2|). 

3.6. Assumption of Local Thermodynamic Equilibrium 

In our derivations we have assumed that the photons and 
electrons are thermally equilibrated throughout the shock pas- 
sage and subsequent expansion. This was done for the sake of 
attaining concrete results and simple analytic expressions, but 
it begs the question of how accurate this is. 

After the initial shock passage, the envelope is heated 
and the electron and photons can reach the same tem- 
perature via Comptonization on a timescale fcomp ~ 
(m e c 2 /kBT(po))/(Kpoc), where m e is the electron mass and 
T(po) is solved from equation ( [T6j >. This timescale must be 
less than the local expansion timescale, f exp ~ Hq/ AV(po), 
otherwise significant expansion occurs before equilibrium is 
reached. For the condition fcomp ^5 t exp , we find that the den- 



sity must be greater than 



po > 3 X 10 (gg/K^Vc 



0.59 T/ 0.030 0.0055 



PI 



gem 



(37) 



This is not much greater than the density of the shock breakout 
(eq. IHl), so we conclude that at least initially all but the very 
outer envelope reaches equilibrium. 

Subsequently, expansion and adiabatic cooling can drive 
the photons and electron back out of thermal equilibrium. 
And in fact, the timescale for thermalization via Comptoniza- 
tion is always shorter than the thermal diffusion time at the 
diffusion depth. Since the photons and electrons were cou- 
pled in the past, we don't expect the spectrum to be altered 
too greatly, but this should be quantified by more detailed cal- 
culations in the future. 
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FIG. 4. — The top panel plots the absolute magnitude of the expanding and 
cooling shock heated envelope. The upper curves show the bolometric lu- 
minosity, and the lower ones show the absolute magnitude in the wavelength 
range of 300—700 nm. The solid lines (dashed lines) are for R» = 3 X 10 s cm 
(6 X 10 s cm). The bottom panel shows the distance out to which an m = 24 
limited optical survey could see such an event. 



4. OPTICAL LIGHT CURVE 

In Figure |4] we show the bolometric and optical (300-700 
nm) absolute magnitude for two different progenitor radii, as- 
suming that the WD emits as a blackbody. In the bottom panel 
we plot the distance out to which such an optical event can be 
seen for an m = 24 limited exposure. The optical peaks at 
around w 1 day following shock breakout. The observabil- 
ity depends sensitively on the WD radius at the time of shock 
breakout. 

Although we have assumed a blackbody emitter for these 
light curve estimates, the atmosphere will in fact be scatter- 
ing dominated. This hardens the spectrum, making the cool- 
ing tail dimmer in the optical. Also, for the wavelengths that 
dominate at « 1 day, metals may play an important role in 
setting the opacity, especially since they will experience some 
recombination for the temperatures at this time. More sophis- 
ticated spectral modeling is needed in order to predict the pre- 
cise shape of the light curve. 

In Figure |4] we have focused on the light curve due to 
from the cooling of the shock-heated envelope, but eventu- 
ally this will be overtaken by the nickel decay. Whether 
or not the shock-heated cooling can be seen above the ris- 
ing light curve of nickel decay depends on the power of the 
early time nickel decay and the WD radius at the time of 
the DDT. If we extrapolate th e typical L oc t 2 law found em- 
pirically (Conl evet al.ll20 06), or the exponential luminosity 
function of lArnetti ([1982), back to early times, they would go 
right through our optical light curves, indicating that nickel 
decay will be comparable to our optical cooling luminosity 
at « 1 day. Since such extrapolations to early times are un- 
certain, early sampling of the SN la light curves is needed 
to better address this issue. As our calculations show, de- 
tecting the signature of a SN la at such early times is pos- 
sible with current and future surveys. Using a SN la rate of 
2.93 x 10" 5 yr" 1 Mpc" 3 dDildav et al.ll2008l) . we estimate that 
w 10-60 SNe la can be observed at these early times each 
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FIG . 5 . — The depth of the thermal diffusion wave during cooling following 
shock-heating of the WD envelope as a function of time. Also marked by dot- 
ted lines are estimates for the maximum mass helium shell 10~ 3 Mq, Iben 
& Tutukov 1989; Shen & Bildsten 2009b) and the maximum mass hydrogen 
shell (« 10~ s Mq; Shen & Bildsten 2009a) that can survive on a 1 .4 Mq 
WD without igniting. This shows what layers of the WD surface may be 
probed by the thermal diffusion wave to help discern between possible SN la 
progenitors. 

year. 

Detailed spectral modeling of the shock-heated cooling 
may also help discern between SN la progenitor models. In 
Figure|5] we plot the thermal diffusion depth in units of mass. 
This shows which mass shells of the progenitor WD are being 
probed by the cooling wave as a function of time. Depending 
on the progenitor model in question, different compositions 
are expected in these surface layers, and once the thermal dif- 
fusion wave is below them, these elements may imprint their 
presence as spectral features. As a comparison, we plot dot- 
ted lines denoting the characteristic maximum thickness of a 
hydrogen shell (» 10" 5 M ; Shen & Bildsten 2009a) or he- 
lium shell (w 10" 3 M Q , Iben & Tutukov 1989; Shen & Bild- 
sten 2009b), as expected to be present in the single degenerate 
scenario for SNe la. 

5. SUMMARY AND CONCLUSIONS 

We have calculated the cooling of a shock heated WD enve- 
lope due to a DDT during a SN la. The general features of the 
resulting flash are a short (~ 10~ 2 s), hard X-ray flash with a 
total energy of ~ 10 40 ergs. The flash is followed by a cooling 
tail as the thermal diffusion wave travels back into the ex- 
panding envelope. The light curve transitions from the X-rays 
to the ultraviolet, and eventually, visual wavebands, reaching 
a peak absolute magnitude of My = -9 to -10 at w 1 day. 
Whether or not this emission is observable as a distinct com- 
ponent in the SNe la optical light curve depends on the WD 
radius at the time of DDT and the nature of the early-time 
nickel decay. Nevertheless, the detection (or lack thereof) of 
emission during the first ~ 20 days before peak would be an 
important constraint on the occurrence of DDTs and on the 
amount of expansion during the deflagration phase. 

A similar shock breakout event may also b e associated with 
the AIC of a WD t o a ne utron star. In fact, iTan et alJ (120011) 
and iDessart et a l. (2007) considered whether the breakout 
shock from core bounce in this case would be a gamma- 
ray burst progenitor, concluding the answer is no. But this 
does not preclude a less powerful, but nevertheless interest- 
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ing, shock breakout as we have described here. The energy 
is greater in the AIC case than the SN la case, because it is 
powered by the gravitational binding energy of the WD (me- 
diated by neutrinos). Using a ty pical energy of 10 50 ergs de- 
posited into 0. 1 Mq (as found by lFryer et al.lll999l) . this gives 
V run w 10 9 cm s" 1 , but with a much larger p mn ks 10 7 g cm" 3 . 
This provides a late time light curve that is brighter by a factor 
of w 2.5, making it easier to observe. Ho wever, the AIC rate 
is much lower than the SNe la rate (~ 1%, Yunge lson & Liviol 



1998), and the shock breakout may be confined along the 
poles (iDessart et al.l f2006. 2007), both of which will make it 
difficult to catch an AIC shock breakout. 



We thank Lars Bildsten, Chris Matzner, and Eliot Quataert 
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